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ABSTRACT 


The  behavior  of  the  system  hazard  rate  function  of  a  two 
component  parallel  system  is  investigated.   The  inter- 
relationships between  the  probabilities  that  the  components 
composing  the  system  are  alive  and  the  system  hazard  rate  is 
examined  with  special  attention  to  certain  points  where 
there  are  important  changes  in  the  behavior  of  the  hazard 
rate  function.   The  behavior  of  the  system  hazard  rate 
function  is  shown  to  depend  upon  the  rates  of  change  of  the 
probabilities  that  the  components  of  the  system  are  alive. 
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I.   INTRODUCTION 

In  the  field  of  reliability,  there  has  been  an  extensive 
study  of  the  properties  of  the  increasing  hazard  rate  (IHR) 
and  decreasing  hazard  rate  (DHR)  classes  of  life  distribu- 
tions.  The  exponential  life  distributions  with  their 
constant  hazard  rates  form  a  dividing  line  between  these  two 
classes.   It  is  reasonable  to  expect  that  the  life  distribu- 
tion of  redundant  systems  composed  of  devices  with  exponen- 
tial life  distributions  fall  Into  the  IHR  class.   However, 
when  two  devices  with  different  exponential  life  distributions 
are  combined  in  a  parallel  redundant  system,  the  system 
hazard  rate  is  at  first  increasing  and  then  decreasing.   The 
underlying  causes  of  this  interesting  behavior  is  the  subject 
of  the  thesis. 

The  usual  definition  of  the  probabilistic  hazard  rate 
function  of  a  device  with  a  random  lifetime  T  will  be  used 
throughout  this  thesis.   The  hazard  rate  r(t)  of  a  device 
is  the  density  function  f(t)  for  the  lifetime  of  the  device 
divided  by  its  survival  function  F(t)=P[T>t]. 

A  distribution  has  an  increasing  hazard  rate  if  r(t)  is 
monotone  non-decreasing  in  t,  and  has  a  decreasing  hazard 
rate  if  r(t)  is  monotone  non-increasing  In  t.   References  1 
and  2  contain  an  extensive  survey  of  the  properties  of  the 
IHR  and  DHR  classes  of  life  distributions.   The  exponential 
life  distributions  are  included  in  both  of  these  classes. 


In  Reference  3»  Esary  and  Proschan  gave  sufficient  conditions 
for  a  system  to  have  an  increasing  hazard  rate  when  it  is 
composed  of  identical  components  with  increasing  hazard 
rates.   A  parallel  system  of  two  components,  with  identical 
exponential  life  distributions,  satisfies  these  conditions 
and  it  has  an  increasing  hazard  rate.   Esary  and  Proschan, 
also  gave  an  example  to  show  that  when  two  components  with 
different  exponential  life  distributions  form  a  parallel 
system,  the  system  survival  function  need  not  have  an 
increasing  hazard  rate.   In  this  case  the  system  hazard  rate 
initially  increases,  overshoots  and  then  decreases  to  an 
asymptotic  value  equal  to  the  lowest  of  the  two  component 
failure  rates. 

In  Section  II,  the  derivation  of  the  system  hazard  rate 
for  a  system  of  two  parallel  components  having  constant 
non-identical  failure  rates  is  reviewed.   In  Section  III, 
a  characterization  of  the  system  hazard  rate  as  a  function 
of  some  conditional  state  probabilities  is  presented  and 
their  interrelationships  are  discussed.   In  Section  IV,  an 
empirical  approximation  for  the  time  at  which  the  system 
hazard  rate  is  a  maximum  is  exhibited. 


II.   SYSTEM  HAZARD  RATE 

Consider  a  system  of  two  independently  functioning 

components  with  exponential  lifetimes  T, ,T2  and  survival 

functions 

-X-,t 
F^t)  =  e   X 

-X~t 
P2(t)  =  e  *         . 

Assume  without  loss  of  generality  that  X?>X-..   The  equations 
for  the  system  survival  function  and  the  system  hazard  rate 
are  generally  known. 

With  the  components  in  parallel,  the  system  lifetime  T 
has  a  survival  function 


F(t)  =  P[max  TlsT2>t] 


-x1t  -x2t  -(x1+x2)t 

=  e    +e    -e 


The  density  of  T  is 


.f(t)  =  -Is-FCt) 

-X,t    -X^t         -(A,+A5)t 
=  Xxe   x  +X2e  d   -(X1+X2)e   L     d 


The  system  hazard  rate  is 

r(t)  =  ti&l 
F(t) 


-x,t    -XQt        -(X  +xP)t 
Xxe   x  +X2e  d   -(X1+X2)e   L      d 

F(t) 


Examining  r(t)  at  t=0  and  as  t  approaches  infinity  gives 
the  expected  results; 

r(0)  =  0 


lim  r(t)  =  min(X, ,Xp)  =  X^ 


The  derivative  of  r(t)  with  respect  to  time  is 


t(t)  =  F(t)f'(t)  +  f2(t) 
F2(t) 

[X^e  x   +  X^e  l   -(X1-X2)£]e   L     d 

F2(t) 


Looking  at  the  sign  of  r'(t),  it  is  clear  that  the 

hazard  rate  is  increasing  if  and  only  if 

2  "*i*    ?  "^p^  2 

X  e    +X   e     >  (X -,-Xp)    and  decreasing  if  and  only  if 

?  "^i^   ?  -Apt  p 

X2  e  L   +  X^e  d      <    (X1-X2)<^  .   If   X^X   ,   r(t)  is 

increasing  for  all  t>0.   The  system  hazard  rate  must  be  at 
a  maximum  for  the  value  of  t  which  is  a  solution  to 

X2  e    +X1  e    -(X1-X2)^  =  0 


III.   SYSTEM  HAZARD  RATE 
IN  TERMS  OF 
STATE  PROBABILITIES 


If  the  parallel  system  is  alive  at  some  time  t>0,  it 
exists  in  one  of  three  possible  states,  i.e.  both  components 
alive,  component  one  alive  and  component  two  dead,  or 
component  one  dead  and  two  alive.   The  system  hazard  rate 
and  its  derivative  have  simple  expressions  in  terms  of  the 
conditional  probabilities  for  these  states  and  their 
derivatives.   The  probability  that  the  system  is  in  state 
i,  1=0,1,2,  given  it  is  alive  at  a  time  t,  will  be  called 
the  state  probability,  denoted  by  P. (t).   For  a  system  with 
lifetime  T  and  component  lives  T_  and  Tp ,  the  state 
probabilities  are: 

(1)  both  components   alive 

P0(t)    =   PCT1>t,T2>t|T>tD 

P[T1>t,T2>t,max(T1,T2>t)] 
"       P[max(T1,T2>t) J 

-(X-j  +  ApH 

e 

F(t) 

(2)  component  1  alive  and  2  dead 

P1(t)  =  P[T1>t,T2<t|T>t] 

P[T1>t ,T2<t ,max(T1,T2>t) ] 


P[max(T1,T2>t)] 


-x,t  -(x1+x2)t 

P  (t)  =  £ =S 

1  P(t) 


-x2t 

=  1  -  e — 


F(t) 


X5t 

=  (e  d   -l)Pn(t)  , 


(3)  component  2  alive  and  1  dead 
P2(t)  =  P[T1<tJT2>t|T>t] 

P[T1<t,T0>t ,max(Tn ,T0>t)] 


P[max(T1,T2>t)] 


-x2t  -(x1+x2)t 

e -e 

F(t) 


-X,t 

1  -  £ 


F(t) 


X,t 
(e  x  -l)PQ(t)  . 


It  is  clear  that  PQ(t)  +  P1(t)  +  P2(t)  =  1,  for  all  t. 

Examining  the  state  probabilities  at  t  =  0  and  as  t 
approaches  infinity  produces  the  following  results. 


P0(0)  =  1,   lim  PQ(t)  =  0 


P, (0)  =  0,   lim  P, (t)  =  1 
P2(0)  =  0,   lim  P2(t)  =  0 
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At  first  glance,  It  might  seem  that  P, (t)  would  equal 
Pfi(t)  at  the  mean  Vrmb  lu  failure  for  component  2,  hut, 
P1(t)  =  PQ(f)  if  and  only  if 

-x  t  -(x1+x2)t     -cx1+x2)t 

e    -e  e 


or 


F(t)  F(t) 


-x1t     -(x1+x2)t 

e      =   2e 


-A9t 

or  1   =   2e  d 


> 


or  t   =  ~  In  2 

x2 


Similarly,  Pp(t)  =  PQ(t)  if  and  only  if 


t  -   7-  In  2 
Al 


The  probability  that  a  particular  single  component  is  alive 
at  some  time  t  is  greater  than  the  probability  that  both  are 
still  alive  for  any  time  t  greater  than  the  mean  tine  to 
failure  of  the  other  component  multiplied  by  a  constant  In  2 

Further,  for  X2>X,,   P-,(t)  >  Pp(t)   for  all  t,  since 
P1(t)  >  Pp(t)  if  and  only  if 

-x2t       -x1t 
l  -  £ >  l  -  S 


F(t)         F(t) 
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or 


or 


-Apt     -At 

e     <  e       , 


O        1     * 


The  derivatives  of  the  state  probabilities  with  respect 
to  time  are : 


p;<t)  - 


-(A,+Ap)t    -(A,+Ap)t   , 
F(t)[-(A1+A2)e    L      d        -e    X      d        F  (t) 


F2(t) 


-(A1+A2)t 

[  "£&1   (x1+x2)  ] 


P(t)        F(t) 


=  PQ(t)[r(t)  -  (A1+A2)]  . 

-A0t     -A_t   , 
F(t)(-Ape   2  )  -  e   2   F  (t) 

Pl(t)  = " =2 

1  F2(t) 


-A2t 

•  [A2  -  r(t)]  . 


F(t) 


-A.t 


1  )  -  e  1      F  (t) 

Pp(t)  = 1 5 

d  F2(t) 


[A,  -  r(t)]  . 

F(t)   1 
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P  (t)  decreases  in  time  from  PQ(0)  =  1  to  lim  Pq^)  =  °>  &nd 

,  t->°° 

Po(t)<0   for  all  t.   P-,(t)  is  increasing  in  time  from 

P1(0)  =  0  to  lim  P1(t)  =  1,  and  p|(t)>0   for  all  t. 
t-*-°° 
Pp(t)  has  a  maximum  at  a  time  x  which  can  be  found  by 

t 

substituting  for  r(t)  in  P2(t)  and  setting  this  expression 

equal  to  zero,  i.e.  when 

-(X1+X2)t      -(2X1+X2)t 
(X-,-X-)e          +  X0e 
__± — £ £ =   o   . 

F2(t) 

(x1+x2)t 

Multiplying  both  sides  by  e         and  combining  terms 
gives 


"~A-i  v  A  ^  ~*  A -. 

e     «- 


and 


X2 


_i_  ln  (__L_) 

X]_      X2~X1 


The  value  of  Pp(t)  at  maximum  is 

X 

1  +  IT 

Aq   —  A-iAq  "~  (j  Aq~A-|  J 

X22  -  X1X2  +  X2X2/X1(X2-X1)  -  (X2-X1)1+(X2/X1) 


and  this  is  less  than  one  for  all  ^2>Xl  * 

As  a  consequence  of  the  definition  of  the  state 
probabilities  and  the  requirement  that 


P0(t)  +  P1(t)  +  P2(t)  =  1 


13 


the  derivatives  are  related  and 


PqU)  +  pj(t)  +  Pg(t)  =  0  . 


Therefore,  Pg(t)  =  -[PQ(t)  +  P-^t)]  and  P2(t)  increases 

t       t  •       t 

when  -?0(t)  >  P-,(t)   and  it  decreases  when  -?0(t)    <   P-,(t)  . 

Pp(t)  increases  as  long  as  PQ(t)  decreases  in  time  faster 

than  P-,(t)  increases,  reaches  its  maximum  value  when 

f       t 
-P0(t)  =  P, (t),  and  then  decreases  because  P, (t)  is  increas> 

ing  in  time  faster  than  Pn(t)  is  decreasing. 

A  hazard  rate  may  be  thought  of  as  a  conditional 

instantaneous  probability  of  failure  at  time  t,  given 

survival  to  time  t.   The  hazard  rate  for  the  system,  given 

that  it  is  in  one  of  its  three  states,  is  0,  A  ,  and  Xp 

respectively.   The  system  hazard  rate  may  be  written  in 

terms  of  these  hazard  rates  and  the  state  probabilities  as 


r(t)  =  0-P0(t)  +  X1P1(t)  +  X2P2(t) 


-(X1+X2)t      -X.^  -(A  +A2)t     -Apt  -(A1+A2)t 

=  0  •  £ +  x1^— f£ +  x2^ =S 


P(t)  P(t)  c    F(t) 

which  reduces  to  the  expression  previously  given  for  r(t). 

The  derivative  of  the  system  hazard  rate  with  respect  to 
time  may  be  written  as 

r'(t)  =  X1PJ(t)  +  A2P2(t)   . 
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Then  the  condition  for  the  maximum  point  of  the  hazard  rate, 
r(t)  =  0,  requires 


X1P|(t)  +  XgPgCt)  =  0   . 


Substituting  for  P1(t)  and  P2(t)  yields 


-X«t  -X,t 

X-, [X5  -  r(t)]  +  \     £ l\     -   r(t)]  =  0   . 

1  F(t)    d  l     P(t)    L 


Solving  this  for  r(t)  gives  its  maximum  value 


rX,t     -Xpt 
X1X2[e   1   +  e  d    ] 

r(t)  =  =0 ^XTt- 

Xxe  d      +  X2e  L 


\ 


e 


-x2t 


From  P,  (t)  = [X0  -  r(t)]  >  0   for  all  t,  it  is  clear 

1      P(t)    d 


that  r(t)  <  Xp   for.  all  t.   At  its  maximum  value   r(t)  >  X- 


if  and  only  if 


or 


-X  t    -Xpt 
X  X  [e   X   +  e   d  ] 

-X  t       -X,t     Al 

Xle      +  X2e 


~A-.  t       -Apt      p  —  A_t       —  A-.t 
XlX2e    +*1*2e     >  Ai  e    +X1X1e   x    , 
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or 


1  2  >  1   ' 


or  Xp  >  X, 


Therefore,  at  its  maximum,   X,  <  r(t)  <  X2 

From  p'(t)  = [A,  -  r(t)]  ,   it  is  clear  that 

2      F(t)     1 

1       ^2 
r(t)  =  A,   at  the  time   t  =  t—  In  (•? ? — )   when  P0(t)  reaches 

J-  A-i        ^o~  ^i 

its  maximum  value.   If  t  >  t  ,  then  r(t)  >  X-,.   r(t)  >  X-.   if 
and  only  if 


x:P1(t)  +  x2P2(t)  >  x1  . 


or  X2P2(t)  >  X1(l  -  P1(t))  , 


or  (X2-X1)P2(t)  >  X1PQ(t)  . 


Since 


X,t 
P9(t)  =  (e  1      -  1)  Pn(t)  , 


then  r(t)  >  X..   if  and  only  if 


At 
(X2-X1)(e  x   -  1)  >  X1      , 

and  for  t  >  t  =  f^-  In  ( ,   »  )   , 
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it  follows  that 

X 


2 


Xnt  If-   ln(X,,-X,  )] 

(X2-X1)(e  X  -1)  >  (X2-X1)(e  l    Al     d   -1   -  1) 


-  <vV(x-^7- x) 


=  Xl 


As  Xp  approaches  X-.  in  value  t  becomes  very  large,  and  if 
X,  =  Xp  =  X  ,   the  system  hazard  rate  r(t)  approaches  X 
asymptotically  from  below.   For   Xp  >>  X-,  ,   t  becomes  very 
small  and  the  system  hazard  rate  is  greater  than  X,  almost 
instantly. 

From  r'(t)  =  X,P  (t)  +  XpPp(t)  ,   r(t)  is  increasing  if 
and  only  if  X,P  (t)  +  XpP  (t)  >  0  ,   and  decreasing  if  and 
only  if   X1P1(t)  +  XpPp(t)  <  0  .   The  behavior  of  the  system 
hazard  rate  is  directly  dependent  upon  the  rates  of  change 
of  the  state  probabilities  P-,(t)  and  P«(t).   r(t)  increases 
as  long  as  P.(t)  and  Pp(t)  increase  and  it  continues  to 
increase  as  long  as  X,P. (t)  >  |XpPp(t)|  .   For 
X1P1(t)  <  |X  Pp(t)|  ,   r(t)  decreases. 

Figure  1  shows  the  state  probabilities  and  the  system 
hazard  rate  plotted  as  a  function  of  time.   As  the  failure 
rate  of  component  two  increases  in  relation  to  the  failure 
rate  of  component  one,  these  curves  are  "pushed  towards  the 
left",  the  system  hazard  rate  increases  rapidly,  and  its 
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maximum  value  Increases.   Pq^)  decreases  to  zero  very 
rapidly  while  P-,(t)  increases  to  one  almost  instantly. 
Pp(t)  becomes  more  peaked  and  it  reaches  its  maximum  value 
sooner. 
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IV.   EMPIRICAL  RELATIONSHIP 
.  BETWEEN  TWO  MAXIMA 

An  empirical  investigation  of  the  relationships  between 
the  state  probabilities  and  the  system  hazard  rate  was  also 
undertaken.   The  system  equations  were  solved  for  various 
values  of  their  parameters  which  permitted  comparison  of 
the  resulting  curves  and  verification  of  the  analytic 
results  of  Section  III. 

The  time  at  which  the  system  hazard  rate  reaches  its 
maximum,  T  ,  can  be  found  by  solving 

,  2  "Xlt    ,  2  "X2t    ,,     ,  ,2        n 
L  e     +  X,  e     -  (A,  -  A-)   =  0 

for  t.   While  this  can  be  solved  numerically,  a  useful 
analytic  solution  for  T  has  not  been  found.   However,  it 
was  discovered  that  the  time  at  which  the  system  hazard 
rate  reaches  its  maximum  is  approximately  related  to  the 
time  t  at  which  P«(t)  reaches  its  maximum  by 


T   =  2t. 
s 


From  the  numerical  results,  it  appears  that  T  is 

X2 

related  to  x  by  some  function  of  a  =  t —  ,  the  imbalance 

Al 
between  the  two  components,  i.e.   T  =  g(a)  •  x.   The  exact 

form  of  g(a)  is  unknown  but  it  appears  to  be  some  function 
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whose  values  are  plotted  in  Figures  2,  3,  and  4.   Note 
that  g(a)  is  never  greater  than  2.17,  nor  less  than  2.0 
for  all  1.001  <_   a  <  100,  and  for  values  of  a  >  15, 
g(a)  <  2.01. 
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